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Abstract 

The transverse momentum Qt distribution is computed for inclusive Higgs boson production at 
■ the energy of the CERN Large Hadron Collider. We focus on the dominant gluon-gluon subprocess 

in perturbative quantum chromodynamics and incorporate contributions from the quark-gluon 
and quark-antiquark channels. Using an impact-parameter 6-space formalism, we include all- 
orders resummation of large logarithms associated with emission of soft gluons. Our resummed 
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, results merge smoothly at large Qt with the fixed-order expectations in perturbative quantum 

Q 

chromodynamics, as they should, with no need for a matching procedure. They show a high degree 
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of stability with respect to variation of parameters associated with the non-perturbative input at 
low Qt- We provide distributions da/dydQT for Higgs boson masses from Mz to 200 GeV. The 
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average transverse momentum at zero rapidity y grows approximately linearly with mass of the 
Higgs boson over the range Mz < mh < 200 GeV, < Qt >— 0.18771^ + 18 GeV. We provide 
analogous results for Z boson production, for which we compute < Qt >— 25 GeV. The harder 
transverse momentum distribution for the Higgs boson arises because there is more soft gluon 
radiation in Higgs boson production than in Z production. 
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I. INTRODUCTION 



The breaking of electroweak symmetry in the standard model (SM) of elementary particle 
interactions is achieved through the Higgs mechanism. In the simplest realization, a complex 
Higgs doublet is introduced, and a single neutral CP-even Higgs boson is predicted. Two 
Higgs doublets are expected at low energies in supersymmetric extensions of the standard 
model. The lightest neutral CP-even Higgs particle h resembles the SM Higgs boson in most 
regions of parameter space of broken supersymmetry Direct experimental searches at 
the CERN Large Electron Positron Collider (LEP) based on the presumption of significant 
decay branching ratio into bottom (6) quarks place the mass of a SM-like Higgs state above 
approximately 115 GeV 2j. An alternative analysis, based on the assumption of Higgs 
boson decay into hadronic jets, without 6-tagging, provides a bound of about 113 GeV |3j. 
A good theoretical description of electroweak data in the framework of the SM requires 
the Higgs boson to be lighter than about 200 GeV [2] . Within the minimal supersymmetric 
standi model (MSSM), the upper bound on the mass of the lightest ffiggs state ts roughly 
135 GeV |5[. The search for the Higgs boson is a central motivation for the experimental 
programs at the Fermilab Tevatron and the CERN Large Hadron Collider (LHC), with 
detection techniques guided by theoretical expectations about its production dynamics and 
decay properties. 

A Higgs boson resembling one anticipated by the SM is expected to be discovered at the 
LHC through various partonic production processes and decays to SM particles jflQ]. These 
include 

• gg -> hX, with h -> 77, h -> W+W, or h -> ZZ ; 

• gg — > tihX, with h — > bb or h — > 77 ; and 

• W + W-(ZZ) -> hX, with h -> W + W~, h -> 77, or h -> t+t~. 

The fully inclusive gluon-gluon fusion subprocess Jl], 0] gg — > hX is the dominant mech- 
anism of production. In this process, production occurs through triangle loops of colored 
(s)particles that couple to the Higgs boson and to gluons. In the SM, the most relevant 
contribution is from a loop of top (t) quarks, with a small contribution from bottom quark 
loops. In supersymmetric theories, the bottom quark loop contribution can be important 
for small uia and large values of tan (3, where uia is the mass of the pseudo-scalar Higgs 



boson and tan (3 is the ratio of Higgs vacuum expectation values. In the MSSM, loops of 
bottom squarks may make an important contribution if the mass of the bottom squark is 
light enough |lo| . 

Semi-inclusive production of the Higgs boson in association with top quarks, gg — > tihX, 
has a relatively low rate because of the large masses in the final state. The semi-inclusive 



weak boson fusion modes [ll| W + W~(ZZ) — > hX require identification of (forward) "trig- 
ger" jets from subprocesses such as qq — ► qqWW — * qqhX. Searches at the Tevatron rely 
on associated production of the Higgs boson with weak bosons, W and Z, and the decay 
mode h — > bb [12]. 

For masses m h of the Higgs boson up to about 150 GeV, the SM Higgs boson produced 
in the gluon fusionprocess should be discovered at the LHC in the rare di-photon decay 
mode h — > 77 jy, [7(. Once data are in hand, the immense 77 background from various 
sources can be determined from the data by a measurement of the di-photon invariant mass 
distribution da/dM^ on both sides of the Higgs resonance. Beforehand, precise theoreti- 
cal calculations of the expected differential cross sections for production of the signal and 



backgrounds [Ha, \i% [15J] are important for quantitative evaluation of the required measure- 



ment accuracies and detector performance. Good estimations of the expected transverse 
transverse momentum distributions can suggest selections in this variable that should im- 
prove background rejection. In this paper, we concentrate on behavior of the Higgs boson 
transverse momentum distribution in the region of small and intermediate values of Qt- 

Above W + W~ threshold but below Z-pair threshold, the decay mode h — > W + W~ has a 
branching ratio approaching 100%. Because there are missing neutrinos in the leptonic decay 
modes of the W's, the invariant mass distribution of the Higgs boson cannot be reconstructed 
directly, but the transverse mass spectrum of the (£ + i~~ jEt) system can be measured. This 
distribution is expected to be broad, and straightforward sideband determinations of the 
background are not possible. However, measurements of the backgrounds in regions in 
which the signal is expected to be absent can then be extrapolated into the expected signal 
region provided that good predictions are available of differential cross sections. 

In the SM, the coupling of gluons to the Higgs boson through the top quark loop is sim- 
plified in the limit of large top quark mass m±. In this limit , the relevant Feynman diagrams 
are obtained from an effective Lagrangian 16j . For inclusive production integrated over 
all transverse momentum, the next-to-leading order contributions in perturbative quantum 
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chromodynamics (QCD) are reported in Ref. |17|. The m t — > oo approximation is valid 
to an accuracy of ~ 5% for rrih < 2m t . Within this approach, the total cross section for 
gg — > /iX, is known to next-to-next-to- leading order accuracy 18 1. 

Computations of the transverse momentum distribution for Higgs boson production at 
leading order in perturbative QCD show that the large m t approximation serves well 
when the Higgs boson mass rrih is less than twice the top quark mass and the Higgs boson 
transverse momentum Qt is less than mt- The next-to- leading order contributions to the 
transverse momentum distribution are computed with numerical integration methods in 

n 

Ref. 20]. The doubly-differential cross section in rapidity and transverse momentum of the 
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Higgs boson is computed analytically at next-to-leading order in Refs. 

When the transverse momentum Qt is comparable to the mass rrih of the Higgs boson, 
there is only one hard momentum scale in the perturbative expansion of the cross section as 
a function of the strong coupling a s , and fixed-order computations in perturbative QCD are 
expected to be applicable. However, in the region Qt <C rrih, where the cross section is great- 
est, the coefficients of the expansion in a s depend functionally on logarithms of the ratio of 
the two quantities, rrih and Qt- A well-studied case is the cross section for gauge boson (W, 
Z, or virtual photon) production in hadronic collisions, at measured transverse momentum. 
The logarithmic term is (a s /vr)ln 2 (Q 2 /Q|) where Q is the mass of the boson and Qt its 
transverse momentum. For Higgs boson production at rrih = 125 GeV, the peak of the trans- 
verse momentum distribution is expected to occur near Qt = 14 GeV at LHC energies, as is 
shown in this paper. Correspondingly, \n 2 (m1/Q T ) ~ 19 and (a s ([i)/Tc) \r?(rri\/ Q T ) ~ 0.7 if 
fi = rrih or ~ 1.1 if fi — Qt- In either case, the relevant expansion parameter in the pertur- 
bative series is close to 1, meaning that higher order contributions in perturbation theory 
are not suppressed. When Qt <C rrih, straightforward perturbation theory is inapplicable. 
Indeed, even at leading order in a s , the computed cross section da/dQ T includes a term 
proportional to (a s /Q T ) \n{rn\/Q T ). This divergence as Qt — > is obviously unphysical. 

All-orders resummation is the established method for mastering the large logarithmic 
coefficients of the ex pan sion in a s and for obtaining well-behaved cross sections at interme- 



diate and small Qt |23l I24L |25| . It is a reorganization of the perturbative expansion that 
brings the large logarithmic terms under systematic control, providing finite predictions for 
the Qt dependence that are different from those at fixed order. Resummation overcomes the 
limitations of fixed-order calculations but still incorporates all the information that these 
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calculations offer. 

Resummation may be carried out in either Qt |25( or impact parameter (b) space, which 
is the Fourier conjugate of Qt space. All else being equal, the b space approach has the 
advantage that transverse momentum conservation is explicit. Using renormalization group 
techniques, Collins, Soper, and Sterman (CSS) [24j devised a b space resummation formalism 
that resums all logarithmic terms as singular as (1/Q T ) ln n (Q 2 /Q T ) when Qt 0. This 
formalism has been used widely for computations of the transverse momentum distributions 
of vector bosons in hadron reactions 26, Higgs bosons [35l 

y,y,y,Q, and other processes [14 l40l]. 

In many of these applications of the CSS formalism, two difficulties significantly limit 
its predictive power in the region of small and intermediate values of Qt- The first of 
these is the ambiguity, if not discontinuity in Qt, associated with the "matching" of the 
resummed and fixed-order calculations in a region of Qt below which the resummed form is 
used, and above which a fixed-order expression is used. The second difficulty is associated 
with the quantitative importance of assumed forms for non-perturbative functions that, 
in many implementations, dominate the Qt distribution in the region of very small Qt 
and affect even the Qt distribution at large Qt- In recent papers 2, Qiu and Zhang 
demonstrate that both of these drawbacks can be overcome. By using an integral form 
for the oscillatory Bessel function in the Fourier transform, they show that the 6-space 
resummation procedure produces smooth distributions for all Qt < Q, obviating the need for 
a matching prescription. Second, to substantially reduce the influence of phenomenological 
non-perturbative functions, they derive a new functional expression that extrapolates the 
perturbatively calculated 6-space resummation expression from the region of small b, where it 
is valid, to the region of large b where non-perturbative input is needed. Their analysis shows 
that perturbation theory itself indicates the form of the required non-perturbative power 
corrections. They also demonstrate that the predictive power of resummation improves 
considerably with total center-of-mass collision energy v^S*- In this paper, we further develop 
the methodology of Qiu and Zhang and use it to study the production of Higgs bosons. 

In Sec. II, we review the pertinent aspects of the 6-space resummation formalism, with 
particular emphasis on Higgs boson production. Working at next-to-leading-logarithm 
(NLL) accuracy, we resum both leading and next-to-leading logarithmic terms associated 
with soft gluon emission to all orders in a s . At LHC energies, the typical values of the 



4 



incident parton momentum fractions xa ~ xb ~ rrih/VS ~ 0.009 (for = 125 GeV) are 
small, and the gluon distribution evolves steeply at small x. Consequently, the saddle point 
in b of the Fourier transform from 6-space to Qt space is well into the region of perturbative 
validity. Resummed perturbation theory is valid at small b, just as perturbative QCD is valid 
at large Qt- However, the Fourier transform from 6-space to Qr-space requires specification 
of the 6-space distribution for all b. In Sec. Ill, we show that the region of large b becomes 
progressively less important quantitatively as y/S increases. The Qt distribution at LHC 
energies is entirely insensitive to the functional form we employ for the non-perturbative 
input at large b, so long as the form used for this extrapolation does not modify the 6-space 
distribution at small b. In Sec. III.B, we discuss the extrapolation from the region of small 
b to the region of large b, and we specify the form we use for the non-perturbative functions 
at large b. We discuss why other approaches in the literature may introduce inappropriate 
dependence on \/S. We provide analytic expressions for the fixed-order perturbative sub- 
processes in Sec. IV. We employ expressions for the parton-level hard-scattering functions 
valid through first-order in a s , including contributions from the glue-glue, quark-glue, and 
quark- ant iquark incident partonic subprocesses. 

Numerical predictions are presented in Sec. V. We show the doubly differential cross 
section da/dydQT as a function of Qt at fixed rapidity y — 0, and the Qt distributions 
for rapidity integrated over the interval \y\ < 2.4. We present results for four choices of 
mass of the Higgs boson, values that span the range of present interest in the SM, from 
m h = M z = 91.187 GeV to m h = 200 GeV. To illustrate interesting differences, we also 
show our results for Z boson production at the same energy. Our resummed results make 
a smooth transition to the fixed-order perturbative results near or just above Qt = Q, for 
all Q, without need of a supplementary matching procedure. Our implementation of the 
resummation formalism works smoothly throughout the Qt region even for the distributions 
differential in rapidity, da/dydQT- Two points are evident in the comparison of Z boson 
and Higgs boson production, with nth = M Z - The peak in the Q T distribution occurs 
at a smaller value of Qt for Z production (4.8 GeV vs. 10 GeV) and the distribution is 
narrower for Z production. The larger QCD color factors produce more gluonic showering 
in the glue-glue scattering subprocess that dominates inclusive Higgs boson production than 
in the fermionic subprocesses relevant for Z production. After all-orders resummation, 
the enhanced showering suppresses the large-6 (small Q T ) region more effectively for Higgs 
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boson production. We compare the predicted Qt distributions for Higgs boson production 
at different masses. The peak of the distribution shifts to greater Qt as rrih grows, in 
approximately linear fashion, and the distribution broadens somewhat. The mean value 
< Qt > grows from about 35 GeV at = Mz to about 54 GeV at = 200 GeV, and 
the root-mean-square grows from about 59 GeV to about 87 GeV. For Z production, we 
find < Qt >= 25 GeV and < Q T > l / 2 = 38 GeV. The harder Q T spectrum suggests that 
the signal to background ratio can be enhanced if Higgs bosons are selected with large Qt- 

Choices of variable parameters are made in obtaining our results, and we examine the 
sensitivity of the results to these choices, including the renormalization/factorization scale 
/i and the non-perturbative input. Scale dependence is the most important source of uncer- 
tainty. It can shift the position of the peak by about 1 GeV, with corresponding changes in 
the normalization of the distribution above and below the position of the peak. The value 
of da/dydQ T at the peak position is shifted by 4 to 5%. Changes in the parameters of the 
non-perturbative input produce effects that at most 1 to 2% depending on the size of the 
power corrections we introduce. In the formulation we use to describe the non-perturbative 
region, there is essentially no effect on the behavior of differential cross section at large 
Qt- In comparison with prior work, we note that the locations of the maxima in the dis- 
tributions da/dydQT occur at slightly larger values of Qt in our case, and the distributions 
themselves differ as a function of Q T above the location of the maximum. The differences 
arise from the different treatment of the non-perturbative input. In our approach, the as- 
sumed parametrization of non-perturbative effects has the desirable property that it does 
not affect the physics in the perturbative region b < 0.5 GeV. 

Conclusions are summarized in Sec. VI. 

II. ALL-ORDERS RESUMMED Q T DISTRIBUTION 

We consider the inclusive hadronic reaction in which a color neutral heavy boson of 
invariant mass Q is produced: 

A(P A ) + B(P B )^C(Q)+X, (1) 

with C = 7*, W ± , Z, or a Higgs boson in the limit in which the top quark mass m t 3> Q/2. 
The square of the total center-of-mass energy of the collision is S. At the LHC, = 14 TeV. 
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In the CSS resummation formalism, the differential cross section is written as the sum 



Arc , (rcsum) aS*) 
U&AB^CX dO AB^CX d<J . 

dQ 2 dydQ 2 T ~ dQ 2 dydQ\ dQ 2 dydQ 2 T 



+ T^lh • (2) 



The all-orders resummed term is a Fourier transform from 6-space 



, (rcsum) 1 r — 

= / tt MQrb) bW A B-*cx(b, Q,x A ,x B ), (3) 

J 271 

where J is a Bessel function. The function WA B ^cx(b,Q,XA,x B ) resums to all orders in 
QCD perturbation theory the singular terms that would otherwise behave as S 2 (Qt) and 
(1/Q T ) hi m (Q 2 /Q T ), for all m > 0. The variables xa and x B are light-cone momentum 
fractions carried by the incident partons from hadrons A and B: 

x A = e y and x B = -0= e~ y , (4) 

and y is the rapidity of the heavy boson. The variables xa and x B do not depend on Qt- 

Resummation treats only the parts of the fixed-order QCD expression that are at least 
as singular as in the limit Qt — > 0. The remainder, including possible less singular 
pieces of the fixed-order perturbative contribution, is defined as the difference of the cross 
section computed at fixed order n in perturbation theory and its Qt <C Q asymptote that 
is at least as singular as Q^ 2 - 

jJX) ^7n- (pert) ^7n- (aSym) 

"°ALl-( \\ (7<T WJ — r.v (lr >MS—( l \ ^ 



dQ 2 dydQ 2 T dQ 2 dydQ\ dQ 2 dydQ 2 T 

This remainder is not significant quantitatively at modest Qt, since the dominant singu- 
larities of the two terms on the right-hand-side cancel in the region Qt — > 0. However, 
the difference becomes important when Q T ~ Q. Explicit expressions for the fixed-order 
remainder terms are presented in Sec. IV. 

We may factor out the lowest order partonic cross section and rewrite the function 
WA B -*cx(b,Q,XA,x B ) that appears in the integrand of Eq. (jHJ) as 

W AB ^cx(b,Q,x A ,x B ) =J2 w ij(b,Q,x A ,x B )alf_^ cx (Q). (6) 
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Here, o'ijLcxiQ) ^ s the lowest order cross section for a pair of partons of flavor i and j to 
produce a heavy boson of invariant mass Q. Because the heavy boson is color neutral, the 
parton flavors ij can be either the quark-antiquark or the gluon-gluon combination. 



A. Higgs boson production 

Specializing to Higgs boson h production in gluon-gluon fusion (and using m h and Q 
interchangeably for the Higgs boson mass in the remainder of this paper), we write 

W A B-+hx{b, Q, x A , x B ) = W gg (b, Q, x A , x B ) ofl_+ hX {Q)- (7) 

In the effective Lagrangian approximation ^| where we keep only the contribution from 
the top quark loop, the lowest order partonic cross section for gg — ► hX is 

^ g l hX (Q)=^^ml6(Q 2 -ml), (8) 

with 

<* - (^) a M < 9 > 

In this expression, Gf is the Fermi constant, and fi r denotes a renormalization scale. We 
note that o"o is formally of second order in the strong coupling strength. However, we define 
o"o be the effective O'th order cross section and focus on the influence of yet higher order 
contributions. If we were to include contributions other than from top quarks in the loop, 
whether from bottom quarks or bottom squarks, we could easily replace <7o by a different 
expression. This inclusion would change the overall normalization of our predictions. For 
on-mass-shell gluons, QCD corrections to the ggh effective vertex are known to order a 2 s 
and likewise could be incorporated in a change in the numerical value of the effective ctq. 

For comparison with Eq. (jHJ), we also present the expression for the lowest order massive 
lepton-pair Drell-Yan cross section: 



<-WxW)=e 9 2 ^- (10) 



r (o) tn\ - .2 
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B. Evolution equation for resummation 



Resummation of the large logarithmic terms is achieved in the CSS formalism through a 
solution of the evolution equation 



d 

;Wgg(b,Q,XA,X B ) 



K g (bpL, a s (ii)) + G g (—, a a (n)) 



Wgg(b,Q,X A ,X B ) , (11) 



<91nQ 2 

where the subscripts g and gg indicate the gluonic process. There are two physical mo- 
mentum variables in this expression, Q and 1/6, as well dependence on the scale variable 
fj, that, in principle, could take on values from of order 1/6 to of order Q. However, the \x 
dependence of K g and G g is determined by renormalization group equations (RGE's), 

d 1 

2^(^,a s W) = —zlgMfi) > ( 12 ) 



<91n/i 



;G g ( — ,a a (fi)) = -~f g (a s (n)) . (13) 



The RGE's for K g and G g ensure that d/d\n(/j, 2 )[Wgg(b,Q : xa, xb)} = perturbatively, an ex- 
pression that is consistent with the fact that W gg is effectively a physical observable. The co- 
efficients 7^ in an expansion of the anomalous dimension 7 s (a s (/i)) = ^„=i lg n Ha: s (fi)/7v) n 
are calculable perturbatively. 

The RGE's for K g and G g may be integrated over ln(/i 2 ) from ln(C 2 /6 2 ) to ln(,u 2 ), and 
from ln(/x 2 ) to ln(C|Q 2 ), respectively, where C\ and C2 are integration constants. The values 
of C\ and C2 are somewhat arbitrary. Following CSS, we select C\ — c — 2e _TE = 0(1), 
where Euler's constant ~ 0.577, and C2 = 1. After the integrations, one derives 

rQ 2 da 2 

Kg(bfi,a s (fi)) + G g (Q/fi,a s (fi)) = - A g (a a (fi)) - B g (a s (Q)) , (14) 

Jc 2 /b 2 [1 

where A 9 is a function of 7 ff (a s (//)) and K g (c, et s (p>)) while 5 g depends on both K g (c, a s (Q)) 
and G g (l,a s (Q)). The choice of Ci = c = 2e~ 7E and C 2 = 1 has the virtue of removing 
logarithmic dependence on these parameters from the functions A g and B g . 

The functions A g and B g are free of logarithmic dependence and have well-behaved per- 
turbative expansions 

Ar(«.(/i)) = E4 B) f^V and ^(« S (^)) = E^ n) f^) n - (15) 

n=l \ 71 / n=l V n J 
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The first two coefficients in the expansions for A and B are known 
choices C\ = c = 2e _TB and C2 = 1, these are 
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For the 



4" = c - 



.(2) _ ^4 



BCD 



n 18 C 
\\C A -2n F 



9 HF 



6 

(l 7T 2 \ 11 „ „ 

C A n F \ 1 C A C F . 

I 12 9 / 8 



(16) 



Here, £(3) is the third Riemann function, n F = 5 is the number of active quark flavors, and 

N 2 -1 4 

C A = N c = 3, and C F = = - (17) 



2N r 



for SU(3) color. The corresponding fermionic coefficients Ay and B«, with i — 1,2 are 



found in Refs. 



HQ 



. The coefficient A^ is available in numerical form 



421. 



The solution of the linear evolution Eq. (jll|) may be written in the form 
W gg (b, Q, x A , x B ) = e~ s ^ W gg (b, - h , x A , x B ). 



(18) 



The function W gg (b, c/b, x A , x B ) depends on only one momentum variable, 1/6, and it can 
be computed in perturbation theory as long as 1/6 is large enough. All large logarithmic 
terms from ln(c 2 /6 2 ) to ln(Q 2 ) are resummed to all orders in a s in the exponential factor, 
with 



S 9 (b,Q) 



rQ 2 dp 2 


In 


(Q 2 \ 


Jc 2 /b 2 p? 




\a 2 ) 



A g (a 3 (n)) + B g (a s (/i)) 



(19) 



The appearance of {c/b) 2 as the lower limit of integration serves to cut off the non- 
perturbative region of small Q F - 

The function W gg (b,c/b,x A ,x B ) may be written in factored form as Q, Q 



c c c 

W gg (b, -, x Al x B ) = f g / A (x A , n, -) f g / B (x B , fi, -) 

where f g / A and f g / B are modified gluon parton distributions 



(20) 



fg/ A (x A ,fi,-) = 



1 d ^ it \r XA c 



(21) 
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Here we see the normal appearance of a factorization scale //, the scale of factorization 
between short-distance hard-scattering and the long-distance physics represented by the 
parton densities. The sum in Eq. (|2T|l runs over all parton types, and 4>u/a(£, A 4 ) is the usual 
momentum (£) distribution of parton a from hadron A at the scale /i. The short- distance 
coefficient functions are expanded in a perturbative series as 

'a s (fi)\ n 



Ca^b (z, (1, ^) = °£lb( Z i b ) 



71 



The first two coefficients in the expansion of C are published (37 



38, 



(22) 



Ci^g(z, /i, 



5(1 -z) 
0. 



C^{z^\) = 6{l-z) 



7T 



°* T + 4 



' \ib* 



-C F z-P^ g (z) In 



(23) 



Here i represents a quark or antiquark flavor, and P g ^ g (z) and Pi^ g (z) are leading order 
parton-to-parton splitting functions, 



■■ 2C A 



1 - ^, . 

l + (l-z) 2 



+ 



11C, 



2n, 



6 



5(1 -z), 



(24) 



There is only one momentum scale, 1/6, involved in the hard-scattering function 
W gg (b, c/b, xa, xb)- Therefore, it is natural to choose \i = c/b in the modified parton densi- 



ties f g (x, /i, c/b) in order to remove the logarithmic terms in Eq. 



However, this choice 



is not required. In principle, /i may be any large momentum scale, e.g., Q. This freedom 
is equivalent to the statement d/d\n(fi 2 ) [W gg (b, c/b, xa, xb)] = order-by-order in powers 
of a s ; or, d/d\rL{i 2 f g (xA,H,c/b) = order-by-order in a s for f g defined in Eq. ()21|) . The 
/^-dependence in C a ^ g is compensated by the //-dependence in the parton densities. To order 
a H , we derive 



d 

dln^i- 



d 

d\n [i 



A(£,aO 



C^l g (-,n, -) 
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+ E 



E/'f 

a Jx S 



d 



d 



,x 



dln^^^^b' 



X . 



M(i - 7 



+ E 

a 

= 0. 



o /a(£,aO 



o:.s 



Pa 



(7) 



(25) 



The last equality follows because both terms on the right-hand-side equal (d/ d In [i 2 )(j)g(x, /i). 

For simplicity, we elect to drop the subscript AB — > hX for the 6-space distribution 
WAB^hx(b, Q, xa, xb)- In its region of validity at small b it takes the factored form 



W*° n (b,Q,x A ,x B )=o 



a,b 



c a . 



4>b/B ® C 6 _ 



x e 



-S{b,Q) 



(26) 



where (g> represents a convolution over parton momentum fractions. Since leading-twist 
perturbative QCD and normal parton densities are valid only for a hard scale \i > fio ~ 1 to 
2 GeV, we expect naively that W pert (b, Q,xa,%b) will be reliable for &<l//iQ~1.0 GeV -1 . 

To perform the Fourier transform in Eq. (|5|L we require an expression for ^(6, Q, x^, ^b) 
that can be used over the entire range of integration in b. To the extent that the region of 
large b is important in the integrand, we must devise a functional form for W(b, Q, xa, %b) 
that extends into the region of large b, a topic addressed quantitatively in Sec. III. B. 

In the CSS formalism, the coefficient functions C^j and the Sudakov form factor S(b, Q) 
are process dependent. As discussed recently in Ref. j^j], it is possible to reorganize the 
perturbatively calculated coefficient functions Cj_»j and the Sudakov form factor S(b, Q) to 
make them universal. This universality is achieved if an all orders pro cess- dependent hard 
part H gg (a s (Q)) is introduced into Eq. (j2H . 



(27) 



W» eit (b, Q, x A , xb) = af^ hX H gg (a s (Q)) £ [<f> a/A <g> CU 9 ] ® [<p b/B ® CU fl ] xe" 5 ^, 

a,b 

with H gg (a s (Q)) = J2 n ^KA)™) and the lowest order if^ = 1. In analogy with parton 
distributions, the universal coefficient functions Cj_>j and the Sudakov form factor S(b, Q) 
can be defined by choosing a "resummation scheme" . The higher order contributions to the 
hard part, can be calculated perturbatively. Since the reorganization affects only the 

perturbative part of W(b, Q, x A , xb), our approach to extrapolate (b, Q, xa, xb) into 

the nonperturbative large b region, discussed in Sec. III. B, could also be used to calculate 
the Qt distribution in this modification of the CSS formalism. 
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III. ROLE OF THE REGION OF LARGE IMPACT PARAMETER 



A. Suppressed role of the region of large b at large y/S 

Prior to presenting alternative expressions for the extrapolation of W(b, Q, x A , Xb) into 
the region of large b, we emphasize one important property of the 6-space resummation 
formalism: the resummed exponential factor exp[— S(b, Q)] suppresses the 6-integrand in 
Eq. Q when b is larger than 1/Q. Therefore, when Q ^> /i , the Fourier transform is 
dominated by a region of b much smaller than l/fio, and the calculated Qt distribution is 



insensitive to the non-perturbative information at 



Saddle-point methods are used in Refs. 



24, 



arge b. 



44] to show that for large enough Q, the 



Fourier transform in Eq. (j3J) is dominated by an impact parameter 



where A = 2A^/(2A^ + (3$) ~ 0.61 and ~ 0.41, for gluonic and fermionic processes, 
respectively, for = 5 quark flavors. This equation shows that the momentum scale 
corresponding to the saddle point, 1/6sp> should be within the perturbative region if Q is 
large enough. In Ref. j^j], a more complete expression is derived for 6 SP that takes into 
account the full b dependence of Wij(b, c/b, x A , %b)- 

The saddle point at Qt = is determined from the differential equation 

-I In (6e" s ^)_ + lln [ 5><^x(Q) W v {b, C -,x A ,x B )) = 0. (29) 



A 



db V Jb=b db 

\v / b=b 

If Wij(b,c/b,XA,XB) has a weak dependence on b near b , the second term in Eq. ()29|) 

can be neglected, and the first term yields bo « bgp. However, as we now show, the b- 

dependence of Wij(b,c/b,XA,XB) is directly proportional to the evolution of the modified 

parton distributions. 

The second term on the left-hand side of Eq. ()29)1 can be expressed as 

d ( c \ d ( c c \ 

— \n{a g g(Q)W g g(b,-,x A ,XB)j = ^ [In f g/A (x A , n, ~) + \nf g / B (xB, V, -)) , (30) 



and 
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d „ ,xa 



d\wb 



E 



d 



<Pa/A((;,V) 
4>g/A{xA,V) ■ 



-P 



71 



a^g\ 



,Xa, 



(31) 



b din fx 2 

The derivative (d/d In fjL)<f>(x, ji) is typically positive (negative) for x < x$ ~ 0.1 (a; > zo), 
and the evolution is very steep when x is far from xq. Correspondingly, the second term in 
Eq. ()29|) should be important when the relevant values of xa ~ Q/vS and 15 ~ Q/\fS are 
much smaller (or much greater) than x . At the LHC, xa ~ %b ~ 0.009 for m fe = 125 GeV. 
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FIG. 1: The function S(b, Q) for Higgs boson production in gluon-gluon fusion at (a) Q = = 125 
GeV and (b) rrih = Mz as a function of impact parameter b. We show results for different orders 
of accuracy in the perturbative expansions of the functions A g and B g . The dotted curve is obtained 
with only the first order contributions Ag 1 ^ and B g 1 ^ retained. For the dashed curve, we include 
A g 2 ^ in addition, and for the solid curve, we keep all terms through Ag and B g 2 ^ in Eq. Mb]) . 

The first term alone in Eq. (J29)) is a decreasing function of b. The derivative vanishes at 
b = b$p producing a first approximation for the saddle point solution, Eq. The negative 
sign in Eq. (|3*T|) , and the fact that the number of small x partons increases when the scale 
\i increases, mean that the second term in Eq. (|2*9*|l is negative when xa and xb are smaller 
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than the typical xo- The full solution to Eq. (|29|) provides a saddle point that is less than the 
first approximation b = &sp- (Conversely, at fixed-target energies where the typical xa and 
xb are greater than xq, the typical b is greater than &sp-) The numerical value for the saddle 
point depends strongly on y/~S, improving the predictive power of 6-space resummation at 
large energies. When Qt > 0, the Bessel function J (z = Qrb) further suppresses the large 
b region of the ^-integration. The suppression is greater for greater Q T , and the 6-space 
resummation formalism is expected to work progressively better for larger Qt- 

In the remainder of this subsection, we present figures that illustrate the b dependences of 
the Sudakov function S(b, Q) and of the factor bW(b, Q) in the integrand of Eq. (JHJ), and we 
comment on their physical implications. The Sudakov function depends only on Q and on 
the perturbatively calculable functions A g and B g , Eqs. (JTHJ-HHJ). Selecting Q = rrih = 125 
GeV, we show the behavior of S(b,rnh) in Fig. ^ The function has a pronounced minimum 
at very small b, b ~ 0.02 GeV -1 . Since it is exp[— S(b, Q)] that enters the integrand, Fig. [T] 
shows that the Sudakov factor strongly suppresses both the large and small regions of b. It is 
instructive to examine the relative importance of next-to-leading order contributions in the 
functions A g and B g . As shown in the figure, it may be sufficient at very small b to retain only 
the first-order contributions, Ag 1 ' and B^\ However, the second-order contributions, both 
A g 2 ^ and are of growing relevance as b grows. The larger Sudakov function at large b in 
next-to-leading order means that the region of small Qt is more suppressed, resulting in a 
shift of the peak of the distribution da/dydQT to somewhat larger Qt at higher orders. The 
large change associated with A g 2 ' may be understood because A g in Eq. (fTT)|) is multiplied 
by \n(Q 2 ). It is possible that inclusion of A g 3 >, once it is known, will produce important 
effects. A comparison of the results in Figs. E^a) and (b) shows that the Sudakov function 
is peaked more narrowly at greater rrih- The suppression of the region of large b is more 
pronounced for larger uih, again influencing a shift of the peak of the distribution da/dydQT 
to somewhat larger Qt as is increased. 

To examine the differences between Z production and Higgs boson production, we show 
in Fig. El the results appropriate for Z boson production, with Q = Mz = 91.187 GeV. The 
second-order fermionic function B^ plays a much less significant role in Z production than 
does the gluonic function B^ in Higgs boson production. Comparison of Figs, ^b) and 
121 shows that the gluonic Sudakov function is much larger than its quark counterpart, but 
the curves are similar in shape otherwise. At b = 1, the ratio is about 2.5. The fact that 
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FIG. 2: The function S(b, Q) for Z boson production in as a function of impact parameter b. We 
show results for different orders of accuracy in the perturbative expansions of the functions A g and 



the gluonic A g and B g functions are larger than their ferrnionic counterparts A q and B q is 
related to the differences in the values of the structure constants Ca and Cp. The larger 
color factor for gluons means that gluons radiate more readily. There will be more gluonic 
showering and correspondingly more suppression of the region of large b (i.e., small Qt) 
when gluonic subprocesses are dominant. Greater suppression of the region of large b means 
that perturbation theory will have greater predictive power and resummed calculations will 
be more reliable for gluonic processes. 

The all-orders resummed 6-space function W(b, Q), Eq. (f2T))). depends on the choice of 
parton densities. We use the CTEQ5M set j^J. In Fig. El we show bW(b, Q) as a function of 
b at the energy of the LHC for two choices of the Higgs boson mass. This function is peaked 
sharply near b ~ 0.05 GeV -1 (c.f., Qt ~ 20 GeV) well within the region of applicability of 
perturbative QCD. The function has essentially no support for b > 0.5 GeV -1 . We expect, 
therefore, that the non-perturbative input at large b will play a negligible role in Higgs 
boson production at LHC energies. Comparison of Figs. H2a) and (b) demonstrates that the 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.61.8 2 

b (1/GeV) 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.61.8 2 

b (1/GeV) 



FIG. 3: The all-orders b-space resummed function bW(b,Q) for Higgs boson production in gluon- 
gluon fusion at yS = 14 TeV for (a) Q = = 125 GeV and (b) = Mz as a function of 
impact parameter b. 

function is peaked somewhat more narrowly at larger m^. 

We present Fig. 0] to illustrate the difference in the resummed b space distributions for 
Higgs boson and Z boson production. In Figs.Efb) andEJ we select the same mass Q = Mz 
and collider energy \fS. The function peaks near b ~ 0.12 GeV -1 for Z production, at about 
twice the value for Higgs boson production. The peak in the Qt distribution will therefore 
occur at a smaller value of Qt for Z production than for Higgs boson production. The 
significantly narrower distribution in b in the Higgs boson case is a reflection of the larger 
gluonic Sudakov factor. Predictions for the Qt distribution of Higgs boson production 
are therefore expected to be less sensitive to non-perturbative physics than those for Z 
production, even if = Mz- (We note that at LHC energies, the typical values of x in 
the parton densities are ~ M z /\/S and ~ m h /\fS. Even for m h = 200 GeV, x < 2 x 10 -2 . 
Therefore, production at the LHC is influenced by the behavior of parton densities at small x, 
not at large x where the more rapid fall-off of the gluon density might otherwise make the Qt 
distribution for Higgs boson production decrease more rapidly than that for Z production.) 
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FIG. 4: The all-orders b-space resummed function bW(b, Q) for Z boson production at = 14 TeV 
as a function of impact parameter b. 

B. Extrapolation into the region of large-6 

The perturbatively resummed function W(b, Q, xa, xb) in Eq. (jSJ) is reliable only when 
b is small. An extrapolation into the region of large b is necessary in order to complete 
the Fourier transform in Eq. (131) . In this section, we summarize different forms used in 

n 

the literature for this extrapolation. The functional form of Qiu and Zhang 32J has the 
desirable property that it separates cleanly the perturbative prediction at small b from non- 
perturbative physics in the large b region. 



In Ref. 

F NP (b,Q,x A ,x B 



24j , a variable 6* is introduced along with a non-perturbative function 



W css (b, Q, x A , x B ) = W peit (h, Q, x A , x B ) F (b, Q, x A , x B ) , 



(32) 



where 6* = b/J 1 + (6/6 max ) 2 < 6 max = 0.5 GeV 1 ; F NP has a Gaussian-like dependence on 
b. The assumed non-perturbative distribution F NP has the functional form [24] 



F i:j (b, Q, x A , x B ) = exp - ln(Q 6 max ) gx{b) - gi/ A (x A , b) - gj/ B (x B , b) 



(33) 
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The ln(Q 2 ) dependence is a prediction of the renormalization group equations. The functions 
9i(b), 9i/A(%A,b), and gj/B(%B,b) are non-perturbative and assumed to vanish as b — ► 0. 
The predictive power of the CSS formalism relies on the derived Q 2 dependence and the 
universality of the functions F?f p , fitted at one energy and used elsewhere. 
Davies and Stirling (DS) proposed a simpler form for F NP 

F% p {b, Q, x A , x B ) = exp [-b 2 ( 9l + g 2 ln(Q6 max /2))] , (34) 

where g\ and g 2 are flavor-independent fitting parameters. To incorporate possible ln(r) de- 
pendence, t = Q 2 / S = xa%b, Ladinsky and Yuan (LY) introduced an additional parameter 
<73 in a modified functional form [3] 

F^(b,Q,x A ,x B ) = exp \ -b 2 {g x + g 2 ln(Qi> max /2)) -bg 1 g 3 \h(100xaXb)] ■ (35) 

Similar to the DS parameterization, no flavor dependence is present. Landry, Brock, Ladin- 
sky, and Yuan (LBLY) reported a global fit to the low energy Drell-Yan data as well as high 
energy W and Z data using both DS and LY parameterizations . 

To preserve the predictive power of perturbative calculations, it is important that 
W css (b, Q, xa, xb) be consistent with the perturbatively calculated W peTt (b, Q, xa, xb) when 
b < &max- However, use of the 6* variable necessarily alters the perturbatively calculated 
6-space distribution even within the perturbative region i32:]. This alteration can be signif- 
icant depending on the choice of the parametrization of F NP . In the last subsection (c.f., 
Figs. El and |U ), we argue that the Qt distributions of W, Z, and Higgs boson production 
at collider energies should not be sensitive to the large b tail of the spectrum. Any signifi- 
cant dependence on the fitting parameters would cast doubt on the predictive power of the 
6-space resummation formalism. 

In Ref. (3^], Qiu and Zhang introduce a new extrapolation 

f W pert (b,Q,x A ,x B ) b<b m&x 
W(b,Q,XA,x B ) = { ' V ' B) wp " max (36) 

[ W^ crt (&max, Q, x A , x B ) F (b, Q; 6 max ) b > 6 max . 

This expression manifestly preserves the QCD resummed 6-space distribution in the pertur- 
bative region at small b. The authors also motivate the adoption of a new non-perturbative 
function in the region of large b: 

F NP = exp { - ln(^p) [ 91 ((b 2 T ~ (&Lx) Q ) + 92 {b 2 - 

-^ 2 (& 2 -6Lx)}- (37) 
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The term proportional to {b 2 ) a , with a < 1/2, corresponds to a direct extrapolation of the 
resummed W pert (b, Q, xa, %b)- The parameters, g\ and a are fixed in relation to H^ pert by 
the requirement that the first and second derivatives of W(b, Q, xa, %b) be continuous at 
b = b mSuX . Dependence on xa and xb is included in the parameters g\ and a. The terms 
proportional to b 2 correspond to inverse-power corrections to the evolution equation. Power 
corrections from soft gluon showers are the origin of the #2 term, and non-vanishing intrinsic 
transverse momenta of the incident partons lead to the §2 term. The values of g% and gi are 
unknown and must be obtained from fits to data. 
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FIG. 5: Ratios of the functions bW(b,Q). Comparisons are shown of (a) the DS parametrization 
and (b) the LY choice with respect to the function we use in this paper. The solid lines show the 
comparisons for Z production at the energy of the Fermilab Tevatron, \/S = 1.8 TeV. The dashed 
and dotted lines show the ratios for Higgs boson production at the LHC. 

In Fig. El we present ratios of the functions bW(b, Q) in use by different groups. In 
Fig. Efa), the numerator is the two-parameter form of Davies and Stirling, whereas in 
Fig. b) it is the three-parameter functional form adopted by Ladinsky and Yuan, Eqs. (p?2*|) 
and However, in both cases we use the updated parameters of Landry et al and the 

CTEQ5M parton densities. The denominator in the ratios is the Qiu and Zhang function we 
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use in this paper, Eqs. (JHUj) an d P7jh with 6 max = 0.5, gi = 0, and #2 = 0. For Z production 
at the energy of the Tevatron, the ratios do not deviate from 1 by more than 20% as long as 
b < 0.5, as might be expected because all groups provide acceptable fits to data. However, 
very significant deviations from 1 are evident at the energy of the LHC, notably in the region 
b < 0.5 where purely perturbative physics should be dominant. The deviations reach 100% 
for b ~ 0.4. This comparison shows explicitly the influence within the purely perturbative 
region of the use of the variable 6* and the nonperturbative function F NP in Eq. (|32p. 

Figure El shows that the deviations from the perturbatively calculated W peTt (b, Q) depend 
strongly on the collision energy y/S. This energy dependence is an artifact of the use of 
the variable 6*. Since 5* < b for all b ^ 0, the shifted 6-space distribution W peTt (b*,Q) < 
W pert (b, Q) for values of b smaller than the location of saddle point, and W pert (6*,Q) > 
W pert (b, Q) for b larger than the saddle point. The magnitude of the difference between 
the W peTt {K, Q) and W pcrt (b, Q) depends on the shape of W pert (b, Q). When yfS increases, 
W pert (b, Q) becomes narrower, and its much steeper behavior leads to a larger deviation 
of W peTt (b*,Q) from W peTt (b, Q). The use of either the DS or the LY parameterizations in 
W (6, Q), Eq. (JH2J), leads to a reduction in the region b < 0.2 and strong enhancement in 
the region b > 0.2. Once the Fourier transform is made to Qt space, the net effects of these 
two differences are a slight shift to smaller Qt in the location of the predicted maximum of 
the distribution do j dydQ? and reduction in the magnitude of the peak value. 



IV. EXPRESSIONS FOR THE FIXED-ORDER CONTRIBUTIONS IN PERTUR- 
BATIVE QCD 

Working at fixed-order in QCD perturbation theory l2lL |22| , one can express the 

differential cross section for Higgs boson production in the form of a convolution 

d°AB^hx _ (o) f 1 f 1 d£ 



i. Jxa C4 Jxa Cfi 



dQ 2 dydQ T "^^Jxa U Jxb Zb 

xH ab ^ hX {QT,Q,^,^-,oi s {fi),fi). (38) 

In this expression, a^^ hX is the cross section given in Eq. (JHj). The functions 4> a /A{xA, 
and 4>b/B{xB, /-0 arc probabilities that incident partons a and b from hadrons A and B 
carry fractional light-cone momenta xa and Xg, respectively, at factorization scale \x. The 
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parton level hard-scattering function H a b^hx(QT, Q, %a/£a, %b/£b, «s(aO, /i) is expanded in 
a power-series in the strong coupling strength a s 



H 



ab~*hX 



/ , n ab^hX 
n=0 



71 



(39) 



The lowest order contribution (n — 0), gg — > h comes from the top-quark loop, and, in 
the absence of partonic intrinsic transverse momentum, it yields a distribution that is a 
5-function in transverse momentum Qt 



Hi ) 



5 1 



? ) 8\Q, 



99 ' K ~ UJ'V & 



(40) 



At first-order in a s , there are virtual loop-corrections to the subprocess gg —> h along 
with contributions in which two partons interact to produce a parton plus the Higgs boson. 
Labeling these 2-to-2 contributions by the tree- subpro cesses , we write the functions H^_^ hX 
in terms of partonic Mandelstam variables s, t, and u as jl9| 
gluon-gluon gg — > hg: 



Hi 1 ) 



<:W 



2tt 



s 4 + t 4 + u 4 + Q i 
siuQ 2 



s + t + u 



(41) 



gluon- quark (antiquark) gq hq (gq — > hq): 



Hi? 



2vr 



C f 



s 2 + u 2 



-tQ 2 



5[s + t + u-Q' 



(42) 



quark (antiquark) -gluon qg — > hq (qg — > hq): 



Hi 1 ) = —C f 



qg 



2tt 



s 2 + g 
-uQ 2 



5(s + i + u-Q' x 



(43) 



and the s-channel quark- antiquark subprocess qq — > /ig: 



(i) 



-(^ fiV 2 
2tt \C\ V 



- 1 



t 2 + u 2 



8(s + i + u-Q' 



(44) 



The first-order contributions from gg, qg, and gq scattering are all singular as Qt — > 0. 
Following convention, we use the word "asymptotic" to designate the terms that are at least 
as singular as Q^ 2 . This asymptotic contribution may be written in analogy to Eq. (J3HJ) as 
the convolution 



da 



(asym) 
AB^hX 



(0) 



dQ 2 algdQ 2 T 



dU f 1 d^ B 



a gg^hX 2—il i I < 

ab j xa C.A J x B SB 



<t>a/A{U, fi) (f>b/B^B,fj) 



XW ab ^ hX {QT,Q, 7^, ^-,a s (fi),fi). 



(45) 
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Expanded in a power series in a s , the asymptotic partonic hard part w a b-^hx becomes 



Wab->hX = 



W 



(n) 

ab^hX 



n=l 



7Y 



(46) 



The lowest order asymptotic contributions (n = 0) vanish for Q T > 0, 



(47) 



Categorizing the first-order contributions as before, we write the asymptotic forms as 
gluon-gluon: 



w 



(i) 

99 



+,5(1 - Z4> Pj-^Zfl) + 5(1 - z B ) Pj-^z^)}, 



where 



, ^b 
-2a = 7—, and z B = — . 
U 4b 



The forms are simpler for the other subprocesses. 
gluon-quark( antiquark): 



w 



(i) 



1 1 



gq 2vr Q\ 



{5(1 - z A ) P q ^ g (z B )} , 



quark( antiquark) -gluon: 



w 



(i) 



1 1 



q9 2ix Q\ 



{5(1 - z B ) P q ^ g (z A )} , 



and quark- antiquark: 



(48) 



(49) 



(50) 



(51) 



w$ = 0. 



(52) 



As discussed in Sec. II, resummation treats only the parts of the fixed-order QCD ex- 
pressions that are at least as singular as Q? 2 in the limit Qt — > 0. The remainder is 



da 



00 



dQ 2 dydQ 2 T 



7 (o) 



1 r 1 d& 



E 7 , , , 

ab ?A Jxb ?B 



V^(^,/i) <t>b/B{£B,H) 



xR a b->hx(QT,Q, 7^, ^-,a s (fi),fi). 
U 4b 



(53) 



As a power-series expansion in a s , the partonic hard part for the remainder term is 
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71=0 \ 

where 

>( n ) _ cr( n ) „„( n ) 



Rab^hX - RatLhX [ S -r~ ) > I 54 ) 



In our numerical work presented in Sec. V, we use order a s expressions for H^]^ hX . In 



n-ab^hX — n ab^hX w ab^hX- \°°J 

ed in Sec. V, we use order a s expressions 
subsequent investigations, one might incorporate the NLO results 0, Q • 

V. PREDICTIONS 

In this section we present our predictions for the transverse momentum distribution of 
Higgs boson production in proton-proton collisions at the LHC energy y/S=U TeV. We 
provide results at both fixed rapidity y = and for the cross section integrated over the 
rapidity range \y\ < 2.4. To illustrate interesting differences, we also show our results for Z 



boson production at the same energy 34]. Our predictions are based on Eqs. (J2J), ©, (J2>J), 
and ([53)1 . We employ expressions for the parton- level hard-scattering functions valid through 
first-order in a s , Eqs. (j4~0J) - (JSJ), with the asymptotic forms Eqs. (}4Tj) - (J52j) . For A 9 and _B g 
in the Sudakov factor, Eq. we take the expansion valid through second order (n = 2), 
Eq. (|lf))l. For the coefficient functions in the modified parton densities, we use the expansion 
through n = 1, Eq. ()23|). (Our results are therefore consistently at next-to- leading- logarithm 
(NLL) accuracy). We use a next-to-leading order form for a s (fi) and next-to- leading order 



normal parton densities jj) |45j. In the fixed-order perturbative expressions that enter 



the "Y" function, we select a fixed renormalization/factorization scale /i = ml + Q^, 
with K = 0.5. For the resummed term, we take as our central value \x = c/b, with c = 2e~ lE . 
To show sensitivity to the choice of fi, we also compute cross sections with the choices 
fj, = 0.5m h and /i = 2m^. 

The extrapolation into the non-perturbative region of large b is accomplished with 
Eqs. ()36j) and (|37|). As remarked before, the parameters g\ and a in Eq. (}3*Tj) are fixed in 
relation to H/ pert by the requirement that the first and second derivatives of W(b, Q, xa, %b) 
be continuous at b = 6 max = 0.5 GeV -1 . Since, as shown in Sec. Ill A, there is little sensi- 
tivity to the region of large b, we adopt for most of our results the values #2 = and ^2 = 
in Eq. (|37j). However, to demonstrate the insensitivity numerically, we also present results 
for other choices of b maiX , g 2 , and g 2 . 
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FIG. 6: Differential cross section da/dydQT for Higgs boson production at y/S =14 TeV and fixed 
rapidity y = with (a) Q = M z , (b) Q = m h = 125 GeV, (c) Q = m h = 150 GeV, and (d) 
Q = rrth = 200 GeV, as a function of transverse momentum Qt- The total prediction is shown as 
a solid line. In (a) and (c), the fixed-order perturbative calculation is shown as a dashed line and 
its asymptotic limit as a dotted line. In (b) and (d), the contribution from the remainder Y term 
is shown as a dashed line. 
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In the limit Qt — > 0, sub-dominant singularities may produce a mismatch [24] between the 
complete fixed-order expressions in Eqs. (j41j) - (}4*4*j) and their asymptotic forms in Eqs. (|48|) - 
(|52|). Sub- leading divergences could drive the difference toward a negative infinity as Qt — * 
0. The "Y" term is not large in the region of small Qt, but the singularity in the limit 
Qt — > can pull the sum in Eq. (J2J) to negative values in the region of very small Qt- There 
are various possible remedies for this difficulty, including introduction of an off-set in Qt 
in the "Y" function. In our numerical work, we use a parameter Q™ m , replacing Qt in the 
"Y" function only by + (Q?! in ) 2 . In keeping with CSS, we select Q^ in = 0.3 GeV as 
our central value and show the sensitivity of results to choices ranging from Q™ m = to 
Qf* = 0.5 GeV. 

The Bessel function in the integrand of Eq. (JHJ) is an oscillatory function, and a high 
degree of numerical accuracy is necessary to ensure accurate cancellations. The number of 
oscillations depends strongly on the value of Qt- When b G (0,2) GeV -1 , the number of 
times Jo(Qt&) crosses zero is 0, 6, 63, and 127 for Qt = 1, 10, 100, and 200 GeV, respectively. 
Following Qiu and Zhang, we use an integral form for the Bessel function 

JJ Z ) = - f cos (z sin(0)) d9 . (56) 

7T JO 

Overall numerical accuracy is controlled by improving the accuracy of the integration in 
Eq. ||5E). 



A. Qt distributions for Higgs boson and Z production 

We display our predicted Qt distributions for Higgs boson production at the LHC in 
Figs. IH1 and respectively. In Fig. we show da/dydQT as a function of Qt at fixed rapidity 
y — 0, whereas in Fig. [7[ we show the Qt distributions for rapidity integrated over the 
interval \y\ < 2.4. (The choice \y\ < 2.4 is motivated by the expected acceptance of the LHC 
detectors for photons We present results for four choices of mass of the Higgs boson, 
values that span the range of present interest in the SM, rrih = Mz = 91.187 GeV, rrih = 125 
GeV, rrih = 150 GeV (where the WW* decay channel is dominant), and = 200 GeV 
(above both WW and ZZ decay thresholds). In all cases, the solid lines represent the total 
prediction, Eq. (J2J). In Figs. Efa) and (c), as well as in Figs. a) and (c), we also show the 
fixed-order purely perturbative result and its Qt —> asymptotic form. In Figs. Efb) and 
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(d), as well as in Figs. Efb) and (d), we show instead the contribution from the fixed-order 
remainder "Y" term. 

In Figs. EJa) and (c), we observe the divergence as Qt — > of the fixed-order purely 
perturbative result and the expected numerical equality of the purely perturbative result 
and its Qt — > asymptotic form at small Qt- The asymptotic form falls away more rapidly 
with Qt- Figures Efb) and (d) demonstrate that the total prediction, Eq. (J2J), is clearly 
dominated by the all-orders resummed term for Qt < Q- At larger Qt, the Y term is 
comparable in size. In the expression for Y, we use the order-a s QCD expressions. Next- 
to-leading order contributions tend to alter theoretical predictions, and, correspondingly, 
we would expect that our predictions would change somewhat for Q T > Q if higher order 
contributions are included. 

The resummed result (labeled "total") makes a smooth transition to the fixed-order per- 
turbative result near or just above Qt = Q, for all Q, without need of a supplementary 
matching procedure. With all Qt distributions extended to Qt = 100 GeV, our imple- 
mentation of the resummation formalism works smoothly throughout the Qt region even 
for the distributions differential in rapidity, da/dydQT- The transition from the resummed 
distribution to the fixed-order perturbative contribution is continuous and less ambiguous 
than is sometimes seen in the literature. 

We note that the resummed results for Higgs boson production begin to fall below the 
fixed-order perturbative expectations when Qt exceeds Q. This effect is likely due to the 
mismatch at order a s between the resummed term and the asymptotic term. A similar 
effect is not seen in our predictions for Z boson production where the fermionic coefficient 
A^p is smaller than the gluonic coefficient Ay\ The fixed-order result should provide the 
correct theoretical answer for Qt ~ Q (and above) where only one large momentum scale is 
present and large logarithmic effects proportional to \n{Q/Q T ) are insignificant. We adopt 
the fixed-order result as our best estimate of the theoretical prediction wherever it is larger 
than the resummed result in the region Qt > Q- 

We display our predicted Qt distributions for Z boson production at the LHC in Figs. El 
The resummed distributions are very smooth functions of Qt over the full range shown 
and they converge without discontinuity with the fixed-order perturbative results above 
Qt = Mz- No supplementary matching is employed. 

Two points are evident in the comparison of Z boson and Higgs boson production, with 



27 




20 40 60 80 100 20 40 60 80 100 

Q T (GeV) Q T (GeV) 



FIG. 7: Differential cross section da/dQx for Higgs boson production at y/S = 14 TeV with rapidity 
integrated over the interval \y\ < 2.4 with (a) Q = Mz, (b) Q = nih = 125 GeV, (c) Q = rrih = 150 
GeV, and (d) Q = nih = 200 GeV, as a function of transverse momentum Qt- The total prediction 
is shown as a solid line. In (a) and (c), the fixed- order perturbative calculation is shown as a dashed 
line and its asymptotic limit as a dotted line. In (b) and (d), the contribution from the remainder 
Y term is shown as a dashed line. 
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rrih = Mz- The peak in the Qt distribution occurs at a smaller value of Qt for Z production. 
At y = 0, the curve peaks at Qt ~ 4.8 GeV for Z production and at Q T ~ 10 GeV for 
Higgs boson production. Second, the distribution is narrower for Z production, falling to 
half its maximum by Qt ~ 16 GeV, whereas the half-maximum for Higgs production is 
not reached until Qt ~ 30 GeV. Technically, both differences are a direct reflection of 
the differences in the functions bW(b, Q) discussed in Sec. Ill A. The 6-space distribution 
is narrower for Higgs boson production and peaks at smaller b. Correspondingly, since 
Q T ~ l/b, the Q T distribution is broader and peaks at larger Q T . The physics behind this 
important difference is that the larger QCD color factors produce more gluonic showering in 
the glue-glue scattering subprocess that dominates inclusive Higgs boson production than 
in the fermionic subprocesses relevant for Z production. After all-orders resummation, the 
enhanced showering suppresses the large-6 (small Qt) region more effectively for Higgs boson 
production. 
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FIG. 8: Differential cross section da/dydQT as a function of transverse momentum Qt for Z 
boson production at = 14 TeV and (a) fixed rapidity y = and (b) \y\ < 2.4. Shown are the 
total prediction as a solid line, the fixed-order perturbative calculation as a dashed line, and the 
asymptotic limit as a dotted line. 
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Comparison of the predicted Qt distributions for Higgs boson production at different 
masses of the Higgs boson shows that the peak of the distribution shifts to greater Qt as 
rrth grows and that the distribution broadens somewhat. At y = 0, the peaks are centered at 
about Q T = 10, 11.5, 12.6, and 14.4 GeV for m h = M z , 125, 150, and 200 GeV, respectively. 
The mass dependence for both Higgs and Z boson production is illustrated in Fig. El For 
Z* boson masses above the physical value Mz we assume that the production model is 
unchanged except for the difference in mass of the Z*. The peak position at y = is at 5.5 
GeV for M z * = 200 GeV. 
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FIG. 9: Location of the peaks of the transverse momentum distributions for Higgs boson and Z* 
boson production as a function of Higgs boson mass at \/S = 14 Te V. For Higgs boson production 
we show results at fixed rapidity y = and for y integrated over the interval \y\ < 2.4. 

The change of the Qt distribution with can also be represented quantitatively with 
plots of the mean value < Qt > and of the root-mean-square < Q\ > 5 > shown in Fig. El 
Because of the long tail of the distribution da /oIQt at large Qt, the mean value < Qt > is 
about 3.6 times the value of Qt at which the peak occurs in the distribution. We observe 
that < Qt > grows from about 35 GeV at = Mz to about 54 GeV at = 200 GeV. 
The curve is nearly a straight line over the range shown, with < Qt >— 0.18rrih + 18 GeV. In 
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computing the averages, we switch from the "total" result for da/dydQT shown, for example, 
in Fig. |U1 to the fixed-order perturbative result when the fixed-order answer is the larger of 
the two choices above Qt = Q- We extend our numerical integrations to the kinematic limit 
in the cases of < Qt > and < Q\ > in order to achieve insensitivity of the results to the 
upper limit of integration. The nearly linear growth of < Qt > and of < Q\ > l l 2 with 
is a reflection of the fact that the peak location in Qt of the distribution da/dydQT grows in 
nearly linear fashion as rrih increases. Since Qt ~ 1/6, the expression for the saddle point, 
Eq. P5|). suggests that < Qt >oc l/6gp oc m\, with fractional power A, rather than a nearly 
linear behavior. We remark that for large rrih, a linear fit is a good approximation to the 
fractional power dependence over a restricted range in m^. 
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FIG. 10: Predictions of (a) the average value of the transverse momentum and (b) the root-mean- 
square for Higgs boson production as a function of Higgs boson mass at \/S = 14 TeV and fixed 
rapidity y = 0. We show results for Higgs boson masses from Mz to 200 GeV. 

For comparison with our results for Higgs boson production, we quote our predictions 
for Z production: < Q T >= 25 GeV and < Q\ > 1/2 = 38 GeV. We note that the difference 
< Q T > — < Qt >— 10 GeV at rrih = Mz, a manifestation of more significant gluonic 
radiation in Higgs boson production. The harder Qt spectrum suggests that the signal to 
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background ratio can be enhanced if Higgs bosons are selected with large Qt, a point to 
which we return in Sec. VI. 

Choices of variable parameters are made in obtaining our results, and it is important to 
examine the sensitivity of the results to these choices. The parameter choices include the 
specification of the renormalization/factorization scale /i, the non-perturbative input, and 
(3™ m . Uncertainties of a different sort are associated with the order in perturbation theory 
at which we work. 
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FIG. 11: Renormalization/factorization scale fx dependence of the differential cross section 
da / dydQx as a function of Qt for Higgs boson mass = 125 GeV at ^/S = 14 TeV and 
fixed rapidity y = 0. We show results for three choices of the scale. In (a) the region of low and 
intermediate Qt is shown in expanded form to illustrate the variation of the position in Qt of the 
maximum of the distribution. 

Dependence on the renormalization/factorization scale /i is often a good indicator of the- 
oretical uncertainty. In Fig. ^2 we show this dependence for the differential cross section 
da/dQT at Higgs boson mass = 125 GeV. Since we are interested principally in exam- 
ining scale dependence of the resummed result, we fix the scale in the purely perturbative 
remainder term Eq. (p)3*|) to be \i — + Q T . This term is not important in the re- 
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gion of modest Qt where the cross section is large. We present results for three choices 
of the scale in the resummed term, Eqs. Q and (|26|). All these results are obtained with 
Q™ n = 0.3 GeV. We select as default choice /i = c/b with c = 2e _TE such that the scale 
varies with the integration variable b. The other two choices are independent of b but are 
proportional to the mass of the Higgs boson, fi = 0.5m^ and fi = 2m^. Scale dependence is 
not insignificant in the region where the distribution is large. It can shift the position of the 
peak by about 1 GeV, with corresponding changes in the normalization of the distribution 
above and below the position of the peak. The value of da/dydQT at the peak position is 
shifted by 4 to 5%. The shift in the location of the peak is about the same as the shift 
that occurs between resummation done at leading-logarithmic (LL) accuracy and at NLL 
accuracy 4f|. 
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FIG. 12: Dependence of the differential cross section da/dydQT on the parameter Q™ in as a func- 
tion of Qt for Higgs boson mass rrih = 125 GeV at ^/S = 14 TeV and fixed rapidity y = 0. We 
show results for three values of Qji m , divided by the results for Q™ m = 0. 

For our central results we select <2™ in = 0.3 GeV. We examine changes associated with 
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this parameter in Fig. The ratio shown is 

da 



R 



shift 



mm- 



da 



(Qr m = 0). 



(57) 



dydQr dydQr 

We vary Q™ m over the range to 500 MeV. The position of the peak of the Qt distribution 
is very stable; within a bin width of 1 GeV, there is no change in the location of the peak. At 
rrih = 125 GeV the change in the value of da/dydQr at the peak location of ~ 11.5 GeV is 
only 0.2% when Q™ m varies from to 0.3 GeV. However, for Qt < 1 GeV, the uncertainties 
are much larger. One message from this exercise is that the resummed Qt distribution 
cannot be computed reliably when Qt <C Qt^, a conclusion that should not be surprising. 
On the other hand, the uncertainties in the Qt distribution in the region of small Qt 
associated with the choice of Q™ m should not be worse than potential uncertainties in the 
data. 
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FIG. 13: Dependence of the differential cross section da j dydQT on the non-perturbative parameters 
(a) £t2 and (b) 6 max os a function of Qt for Higgs boson mass = 125 GeV atVS =14 TeV 
and fixed rapidity y = 0. We show results for two choices of each parameter. 

In Fig. EI] we examine the dependence of our predictions on the choice of the non- 
perturbative parameters & max , 92, and g 2 - We combine the dependence on the latter two 
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parameters in a function G2 = g 2 \n(Q 2 b 2 nSLX / c 2 ) + g%, denned in Ref. |32j. In Fig. ITBTa) we 
plot the ratio of the Qt distribution with residual power corrections over its value without 
these contributions. In Ref. Q, G2 = 0.4 GeV* is shown to produce the best fit to Tevatrou 
data on W and Z boson production. To enhance the possible influence of power corrections, 
we choose parameters that are a few times larger than are needed in W and Z production. 
In Fig. ITBT a). the values G2 = 0.8 GeV 2 and 1.6 GeV 2 are used, twice- and four-times the 
values determined earlier. The corresponding changes in the predictions for Higgs boson 
production are very small. In Fig. ITBT b). we examine sensitivities to the choice of 6 max . The 
ratio displayed is R(b maiX ) = da / dydQT(b mSLX ) / da / dydQT(b ma , x = 0.5 GeV -1 ). For changes 
in frmax between 0.3 and 0.8 GeV -1 , the uncertainties in our predictions are much much less 
than 1%. They are consistent with the level of accuracy of our numerical computations. 
Evident in Fig. is that the formulation we use to describe the non-perturbative region 
has no effect on the behavior of differential cross section at large Qt- 

Comparison of Figs. E] and EH shows that uncertainties related to scale variation are 
greater than uncertainties associated with non-perturbative physics. Uncertainties associ- 
ated with physics in the region of large b are at most 1 to 2% depending on the size of the 
power corrections we introduce, whereas scale dependence can change distributions by more 
than 10%, if one examines the spectrum as a whole, and by about 5% in the region of the 
peak. Non-perturbative effects vanish very quickly as Qt increases, while sensitivity to the 
scale choice increases as Qt increases. Variations associated with non-perturbative physics 
in the large b region are much smaller at LHC energies than higher-order corrections to the 
perturbative functions A, B, and C. 



A different resummation scheme, such as discussed in Ref. 1431. can lead to uncertainties 
similar to those shown in Fig. ^2 of this paper and in Fig. 7 of Ref. [46]. Because our method 
for extrapolation into the nonperturbative large b region is different, we expect the locations 
of the maxima in the distributions da/dydQT in Fig. 7 of Ref. |46j to shift to somewhat 
larger values of Qt- 

VI. CONCLUSIONS 

Discovery of the Higgs particle is essential to shed light on the mechanism of electroweak 
symmetry breaking. The partonic subprocess g+g — > hX dominates Higgs boson production 
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in hadronic reactions when the Higgs boson mass is in the expected range rrih < 200 GeV. 
The two-scale nature of the production dynamics, with mass rrih and transverse momentum 
Qt both potentially large, and the fact that the fixed-order perturbative QCD contributions 
are singular as Qt — > 0, necessitates all-orders resummation of large logarithmic contribu- 
tions in order to obtain meaningful predictions for the Qt distribution particularly in the 
region of modest Qt where the cross section is largest. We perform this resummation of mul- 
tiple soft-gluon emissions using an impact parameter 6-space formalism. At LHC energies, 
the typical values of the incident parton momentum fractions xa ~ %b ~ rrih/yS ~ 0.009 
(for rrih = 125 GeV) are small, and the gluon distribution evolves steeply at small x. Conse- 
quently, the saddle point in b of the Fourier transform from 6-space to Qt space is well into 
the region of perturbative validity. The resummed Qt distributions are therefore determined 
primarily by the perturbatively calculated 6-space distributions at small b, with negligible 
contributions from the non-perturbative region of large b. Our resummation formalism 
has excellent predictive power for the Qt distributions of Higgs boson and Z production 
at the LHC energy, v^S* = 14 TeV, for Qt as large as the respective boson masses Q. 
An examination of the numerical sensitivity of our results to variations in the renormal- 
ization/factorization scale fi and to parameters associated with the form of the assumed 
non-perturbative function at large b shows the expected modest dependence of \i and essen- 
tially no dependence on the non-perturbative parameters, except in the region of extremely 
small Qt- 

In this paper, we present predictions for the Qt distributions of Higgs boson and Z 
production at y/S = 14 TeV. Results are shown in Figs. El-El At the same mass, rrih = Mz, 
the predicted mean value < Qt > is about 10 GeV larger for Higgs boson production than for 
Z boson production. For the Higgs boson, < Qt > grows from about 35 GeV at rrih = Mz 
to about 54 GeV at rrih = 200 GeV, and the root-mean-square < Q T > 1//2 from about 
59 GeV to about 87 GeV. 

Searches for the Higgs boson in its decay to two photons require an excellent under- 
standing of the production characteristics of both the signal and backgrounds. The primary 
background to direct photon production comes from decays of hadrons such as 7r°, 77 , and 
uj°, themselves produced copiously at large transverse momentum from the fragmentation 
of jets, as well as from misidentified photons (e.g., from processes that produce two jets or 
a jet and a photon and in which jets are misidentified as photons). Even after excellent 
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photon/jet discrimination and application of photon isolation restrictions, simulations show 
ratios of the expected signal to irreducible backgrounds in the range of 4 to 5 xlO -2 for 
rrih in the interval 120 to 130 GeV Q. The irreducible backgrounds arise from the QCD 
annihilation qq —>■ 77, the Compton qg —>■ 775 (with one of the final photons produced from 
parton fragmentation (generalized bremsstrahlung)), and the gluon box on — > 77 subpro- 
cesses, and their counterparts at higher orders of perturbation theory [3, [lj, [l^ . Although 
there is some dependence on the di-photon invariant mass and photon transverse momen- 
tum, the contributions of the annihilation and gluon box processes are comparable, and the 
fragmentation processes make a smaller contribution than either of these others. Of interest 
to us is the expected dependence of the ratio of signal to irreducible background. In this 
paper we provide predictions for the transverse momentum distribution of the signal. The 
gluon box subprocess has the same initial state as the inclusive gluon fusion process that 
produces the Higgs boson. Consequently, for di-photon invariant mass M 77 = m^, we expect 
that the transverse momentum spectrum of the gluon box contribution to the irreducible 
background will have the same shape after soft gluon resummation as that for the Higgs 
boson. However, the other major component of the irreducible background, the annihilation 
qq — > 77 subprocess, has the same initial state structure as Z boson production. As we 
show in this paper, the Qt spectrum of Z production is predicted to be softer than that 
for Higgs boson production (c.f., Figs. I§1 and ITU]) . We suggest therefore that a selection of 
events with large Qj 1 could help in significantly improving the signal to background ratio. 
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